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Abstract 



We investigate the effective potential for a scalar theory with spon- 
taneous symmetry breaking at finite temperature. All 'daisy' and 'super 
daisy' diagrams are summed up and the properties of the corresponding 
• gap eqation are investigated. It is shown exactly that the phase transition 

is first order. 



>< 

^ : 1 Introduction 

The investigation of quantum fields at finite temperature is an important prob- 
lem for modern cosmology and particle physics. Since the pioneering papers by 
Kirzhnits f^] and Linde numerous investigations devoted to the development 
of the formalism and to various applications appeared. The best known papers 
are the famous ones by Weinberg and Dolan and Jackiw [Q. There are ex- 
cellent books i, I, 0, i and reviews [|, 0, |n|, |l|, |l|, |14|, [l| on these topics. 



However, the complete understanding of the spontaneous symmetry breaking in 
the early universe remains a key problem. It becomes more important as now 
all particles of the Standard Model (except for the Higgs) are discovered and 
their characteristics like masses are known. Hence, a quantitative analysis of the 
electroweak phase transition can be made which allows a direct comparison with 
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the cosmological constraints among which the most important is the possibihty 
of the baryogenesis scenario |jl5[ , [0 . There is a lot of papers devoted to these 
problems (see, for instance, [17], |18|, and references therein), where the 



temperature effects are taken into account with different accuracy. Of course, in 
order to make such calculations reliable, the theoretical foundations of the for- 
malism must be deeply understood. But there are some subtle problems to solve 
the present paper is devoted to. 

The main problem in field theory at finite temperature with symmetry break- 



ing is the presence of long range correlations resp. infrared divergencies [§, |T2|, [10 



Tj]. To some extent, the perturbative approach itself became questionable. The 
point is the necessity to go beyond perturbation theory, i.e., beyond the simple 
expansion in powers of the coupling constant. One has to sum up certain classes 
of diagrams and to show that thereafter no infrared (or other) divergencies are 
left. There is a large number of papers trying to accomplish this task. Quite 
early it was understood that the so called 'daisy' diagrams must be summed up. 
Already in |Q it was said that even more is necessary - one has to consider the 
corresponding gap equation (an approximated version of the Dyson equation). 
However, to our knowledge, up to now no complete investigation of this kind had 
been done. For instance, the problem with the imaginary part of the effective 
action was not solved completely. Frequently, only the high temperature approx- 
imations of the corresponding quantities had been used for that investigations 
(see for example [T^, |18[|). Although the order of magnitude is captured right by 
this approximation, there is a reason to expect that a more correct calculation 
at finite (instead of high) temperature may change the results for about 10 ... 
20 % which may be important in cosmological estimates. The gap equation was 
considered in including a higher self energy graph in the gap equation. A 
first order transition was found, however with unwanted imaginary parts. Similar 
results have been found in |^ . 

There is another approach using flow equations based on the averaged action 
method (||22| and subsequent papers) where the phase transition was found to be 
of second order. Concerning the order of the phase transition the same result was 
obtained in |2^ using standard renormalization group equation in 0^3^ and 0(^4) 
theories. However, in both approaches there is no sufficiently strong estimate for 
the neglected contributions. 

In the present paper, in order to investigate these problems, we consider the 
scalar $^ theory in 4(=3+l) dimensions with spontaneous symmetry breaking 
at finite temperature in Matsubara representation. Using functional methods we 
sum up all 'daisy' and 'super daisy' diagrams and all 'tadpole diagrams' as well. 
For the effective mass we obtain the corresponding gap equation. All important 
properties of its solutions are derived exactly, approximative methods and numer- 
ical calculations are given. The effective potential is shown to have a first order 
phase transition, all characteristic quantities are calculated, numerically where 
necessary. The imaginary parts of the effective mass and the effective potential 



2 



are shown to be absent in the minima of the effective potential at any tempera- 
ture. So the summation made turns out to dehver the essential properties of the 
effective potential correctly. Further corrections may add powers in the coupling 
only. The methods used are in fact transparent and quite elementary, however 
some background in functional methods is helpful. 

To the authors knowledge, this result had not been obtained before. All 
previous works have been done in one or other weaker approximation, frequently 
using high temperature expansions. We discuss the most important of these 
approximations in detail showing their relations to the correct values. As it 
frequently happens, the correct approximation is more elegant and technically 
simpler that the intermediate ones. 

In the next section the necessary functional formalism is introduced and the 
summation of the infrared dangerous graphs is made. In the third section which 
constitutes the main part of the paper the gap equation is investigated and all 
important properties of its solution are derived. In the fourth section the effective 
potential is investigated in detail and the phase transition is shown to be of first 
order. In the following section the comparison with different approaches is made. 
The results are discussed in the conclusions and some basic functions used are 
collected in the appendix. 



2 Perturbative field theory in functional formu- 
lation 

The appropriate formulation of perturbative quantum field theory is the func- 



tional formulation. We follow the book EM which gives the most complete treat- 



ment of this formalism. Similar representations are given in 0] and in |25| 



Let us first consider a general Euclidean scalar theory with the action 

^ ' 2 3 4 ^ ' 

We use condensed notations. The arguments of all quantities like the field $ are 
dropped. Wherever arguments appear they are integrated over. The first term 
in the action is the 'free part', quadratic in $. Its kernel is the operation K. In 
our case it is 

K = n-m^ . (2) 
In complete notation the action reads 

S{^) = Jdx (^^^x){a - m2)$(x) - - . (3) 

The generating functional of the Green functions is given by 

Z{J) =N Jd(^ exp - - + $7^ , (4) 
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where is a normalization constant and J {J{x)) is a source. We consider it as 
a functional argument. The perturbative representation of Z{J) is 



Z{J) 





'15^5' 




exp 


2 5$ 5$ 


exp 



7.^3 



A 



3 4 



(5) 



where ^ denotes the variational derivative. By expanding the exponentials the 
last formula turns into the representation as the sum over all graphs with lines 
A = —K~^ and generating vertices — f '^''^ — j^'^ + ■ An actual vertex factor is 
the corresponding functional derivative of the generating vertex. A line is 'hung 
up' by the variational derivatives between two (or one and the same) vertices. 
The generating functional of the connected Green functions is known to be 

W{J) = \yiZ{J). (6) 

To pass to the effective action r(<l>) one has to perform the Legendre transforma- 
tion. For this reason one introduces the new functional argument $ = jjW{J), 
expresses J in terms of $ and defines the effective action by 

r($) = W{J) - M. 

Its perturbative expansion is known to be 



(7) 



1. 



r($) = + ^Tr In A + W^^\j A^^^) 



(8) 



Here, 5* is the action (|ip, the 'Tr In' is the one loop contribution and W^^^ 
denotes the one particle irreducible (IPI) graphs contributing to W, by means of 
J A~^$ amputated for their external legs. 

Now the effective potential is minus the effective action for zero argument^]: 



eff 



-r(o) . 



(9) 



In the graphical representation it contains the tree contribution S{0) (which is 
zero for the action (|l|)), the 'Tr In' as the one loop contribution and ly^^^(O), the 
IPI vacuum graphs of W. 

In general, in the representation by graphs, lines closed over one vertex, C') , 
do appear. In field theory without external fields or without symmetry breaking 
resp. at T = these result in unessential constants and are usually hidden in a 
normalization. However, for finite temperature they deliver nontrivial contribu- 
tions and must be taken seriously. In fact, from the considered <l>^- vertex they 

appear with two external legs, >( . These insertions are called daisy graphs. 

^Usually this notation is preserved for a constant background field. Up to now the formulas 
are completely general and the background is not specified. 
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They are like mass insertions. As such they cause the known infrared problems. 
Consequently, they must (and, in fact, can) be summed up completely. The same 
holds for tadpole insertions, ^ , appearing from vertices. 

The problem of summation had been tackled repeatedly. Difficulties arise 
with respect to the correct count of the symmetry coefficients and beside. Here 
we use a simple functional method to sum up these insertions completely. From 
the point of view of the functional methods as described in the book this is an 
easy exercise consisting in adding and subtracting three terms in the action, i.e., 
in redistributing contributions between the 'free part' and the remaining part, 
which is treated perturbatively. So we rewrite the action S (|1|) in the form 



S 



1*( 

9 V 



□ 



m 



(10) 



|a^<l>^ + 6^ and M3 = + c$ are 



where the new vertex factors Ai^ = 
introduced. For the moment the parameters a, b and c are arbitrary constants. 

In the functional language all lines closed over one vertex result from the 
operation of 'hanging up lines' acting only on that vertex. So, for instance, the 
daisies on one $^-vertex are generated by 



exp 



i— A— 

25$ 5$ 



$4 = $4 + 6Ao$^ 



3A^ 



(11) 



where Aq is A for coinciding arguments (see (O)). A detailed explanation can be 
found in p4[. Nevertheless, these formulae call for a more detailed explanation 



here. First of all, the line A has two ends and, hence, it is the function of two 
arguments. For example, in free space it is nothing else than the propagator 
A(x,y) = J The operation with the functional derivatives reads 



- — A — 

25$ 5$ 



da; / dy 



-A{x,y)- 



The vertex is $^ 



5$(x) ' '5$(y) ■ 
/ dx $(x)'^ and the closed line is represented by 

dA; 1 



Ao = A{x,y) 



(27r)4 P + m2 



(12) 



(for the moment we ignore the convergence problems). The index '0' is used 
to symbolize that the line A{x,y) is taken at x — y = 0. Now, in the explicit 
notations, the derivatives can be carried out trivially and returning back to the 
condensed notations the rhs. of Eq. (|l^) will be obtained. Analogously, the 
relation 

^ ^ ■ ^ <^3 + 3Ao<l> (13) 



exp 



A — 

2 5$ 5$ 



can be checked. It describes the emergence of the tadpole diagrams from a 
vertex. 
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In order to avoid the appearance of 'daisies' the constants a and b introduced 
in Eq. ([T0|) may be adjusted accordingly. In demanding 



exp 



25$ 5$ 



Ma = — 



(14) 



this goal can be achieved. At this point it must be noted that by 'adding zero' 
in Eq. (|TOD we changed the first part of the action which is the 'free part' giving 
rise to the lines (compare Eqs. (^ and (|D). Hence the mass in the propagator is 
changed according to m? + . We denote the line with this changed mass 

by a tilde. 

From the formula (|TT]) together with the simpler one 



exp 



1— A— 

25$ 5$ 



$2 ^ $2 



An 



we obtain 



(15) 



exp 



2 <5$ ^ (5* 



exp(A<z 



exp 

-,2 



2 (5<I>"(5<I> 



= + (-|AAo + 4) $2 _ |AA2 + f Ao + 62 
In order to fulfill the requirement (|I4|) the equations 

= 3AAn 

and 



^0 



= --AA 



must be satisfied. 



In a similar way we obtain from (13) 



exp 



2 5* " (5<I> 



exp(A^3) = exp 



2 5$ ^ (S* 



and from the requirement 



exp 



the parameter c follows to be 



25$ 5$ 



3 A 

c = -7A0 . 



(16) 
(17) 



(19) 



So, by virtue of (|Tj) and ([T8|) , no closed lines will appear from the vertices 
and A^4. This is a quite strong consequence because together with the 'daisy' 
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diagrams also all 'super daisy'0 diagrams appear to be summed up. The same 
holds for the 'tadpole' diagrams. This procedure is similar to the concept of the 
'reduced vertex' in p. 38. 

The meaning of Eq. (|16]) is well known - it is a special case of the Dyson 
equation, also called 'gap equation'. We shall follow the latter terminology. To 
illustrate this we denote the mass of the line as explicit argument of the function 
Aq as Ao(M), where = + a^. In this notation the A without tilde is 
Ao(m). So we drop the 'tilde' in what follows. Than Eq. ([TBD becomes 



= m^ + 3AAo(M) . (20) 

M has the meaning of the effective mass. 

We note the graphical representation of Ao(M), 

/dA^ 1 

Now we turn to the corresponding theory with temperature using the Mat- 
subara formalism. Thereby the formulae given above remain unchanged. Further 
we introduce the spontaneous symmetry breaking by giving the mass term the 
'wrong sign'. Than, after m? —m? the field $ must be shifted 

$(x) = + 77(0;) , (22) 

in order to quantize in the vicinity of the lower lying minimum of the action, 
where is the condensate field, a constant, and ri{x) is to be quantized. After 
this substitution the action reads 

S = ^<l>e' - ^<l>/ + - \<^,')v + lv(^-^,')r^- A<l>er/3 - ^r/^ , (23) 



where /i^ = — + 3A$c^ is the new mass parameter squared. Using this action 
instead of (|l|) and performing the summation of the 'daisy' graphs as explained 
just above the action takes the form 

S{v) = ^<fe'-^$e' + ^AAo(M)2 + $,(-m2 + A<l>,2 + 3AAo(M))7^ 

+ir/ (□ -M')r]- A<I>,A^3 - ^M^ , (24) 

= r?^ + 3AAo(M)r/ (25) 



where 
and 



M4 = + G\Ao{M)r]^ - ^XAo{Mf (26) 



^They are also called 'foam' diagrams in EC 
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are the vertices which do not obtain closed loops over themselves in the pertur- 
bative expansion assuming the equation (|20|) holds which reads 

= -m^ + 3A<I>,2 + 3AAo(M) (27) 

in the present notations. 

According to the general formula (11), the effective action takes now the form 



T{v) = -V,s + ^c{-m^ + X^c^ + 5XAo{M))7] + ^7]{a-M^ 



V 



where 



-A$,A^3 - ^M, + W'^\ri) , (28) 

Kff = -^<l>c' + ^$c' - IxAliM) + V1(M) - W'''\0) (29) 
is the effective potential and 

Vi = -- Tr In A 
2 

is the one loop contribution with the mass M in the propagator, see Eq. (|53|). 

In the effective action r(?7), (p8|), as well as in the action (p4D, the term linear 
in r] acquired an additional contribution. The effective mass M in the addendum 
quadratic in 77 must obey the gap equation (^) and becomes a function of the 
condensate $c and of the temperature in this way: M{^c] A, T). All higher graphs 
with external legs are contained in VF^^^. 

In the effective potential Ves the first two terms are the so called 'tree' ap- 
proximation. The corresponding minimum of the effective potential is reached at 
$c = '^Z \/A which is at once the condition for the linear term in T to vanish at 
tree level and the mass of the field is = 2m^, then. 

The tM.d t„o, in can be represented graphically as g looking like a 

'eight'. The function Vi is a notation for the 'Tr In' contribution and ly^^^(O) are 
the higher loop graphs without external legs (vacuum graphs). We do not need 
to consider them in this paper. 

The effective potential has to be considered together with the gap equation 
(pTl). This allows to exclude Ao(M) from explicitely appearing in the expression 

K,._(^£l±i!!!)!,^*._^*.,V,(M), ,30) 

which, together with (0), will be used in the following. 

The functions Aq{M) and Vi(M) can be written down quite explicitely. This 
had been done in a number of papers. Different representations are possible. We 
collect the two most useful ones in the Appendix together with their expansions 
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for high T which is at once the expansion for small M. These quantities con- 
tain ultraviolet divergencies which are removed by a standard procedure which 
includes also the fixing of the normalization conditions. They are chosen not 
to change the minimum of the effective potential known at tree level in case of 
zero temperature (the finite temperature contributions are ultraviolet finite and 
unique). So we require 

^— = and — — ^ = . 

In fact, the superficial divergence degree of vacuum graphs is four. Because Kfr 
depends on rather than on $c, there are only three independent normalization 
conditions to impose. After the first two there remains one condition fixating the 
value of Kfj at some value of <I>c^- It is natural to demand Kff(^c = 0) = 0. 
But the calculation of the corresponding constant to be subtracted from V^g is 
not easy. So we leave that condition open. In fact, it becomes important only in 
figure ^ where it had been adjusted numerically. 

The single 'daisy' graph Aq, which is in fact the first graph in the self energy, 
has superficial degree two and two normalization conditions have to be imposed. 
In the same aim as above we choose 

A(M2)|*.=Wv^ = 0. 

The second condition results in an additive constant. We choose it so that the 
equation 

^^''''^''-A„(M) (31) 



M DM 

holds. Note that this relation holds for the temperature dependent part auto- 
matically. 



3 Solving the Gap Equation 

The solution of the gap equation is the key element of the whole procedure. There 
is, of course, no explicit solution. However, the most important properties can 
be derived exactly. 

In solving the gap equation (^) we obtain the effective mass M appearing 
after summing up the 'daisy' diagrams. It is a function of the condensate $c, 
further it depends on the bare mass m, on the coupling A and on the temperature 
T: M($c; A, T). The inverse function, (^c{M'ifT^-,^)T) is given explicitely by 
Eq. (^: = +w? - 3AA(M))/3A. This function can be easily plotted, 
see Figure |I|, for real values of M. As seen, it is not monotone and the inverse 
function, i.e., M($c; A, T) is not unique. So the gap equation has two solutions 
in M for each $c- In order to understand that better we consider the graphic 
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20 ■ 



15 



Figure 1: The function $c^(Af^) as given by the gap equation (|27|) for A = 5 and 
T = 1 

solution of Eq. (pTl). For this we write the equation in the form = rhs{M) 
with the notation rhs{M) = -rn? + 3A$c^ + 3AAo(M) and plot the Ihs. and the 
rhs. of this equation as functions of in one picture, see figure |^. The curve 
for the rhs. starts at M = in rhs\f^J^^■^ = —m?5i + 3A$c^ + AT^/4 which follows 
from (^) where we introduced the notation 

The derivative is negative there (see the second term in the rhs. of Eq. (0)) so 
that for -m^5i + 3A<l>c^ + ATV4 > there is a solution close to M = 0. In the 
opposite case there is no real solution there. For large values of M, the curve 
describing the rhs. will be dominated by the logarithm, see formula (0). For 
3A ln(M^/2m^)/167r^ 1 it becomes larger than the Ihs. and delivers the other 
solution, seen in Figure 1 as the decreasing part of the curve. It takes place for 
large values of A and M (at any T) and is of no relevance for physical applications. 

For —m?5i + 3A<l>c^ + AT^/4 = the curve of the rhs. starts in the origin and 
the solution is exactly M = 0. 

In general, the solution of the gap equation is non perturbative and can be 
obtained only numerically. However, for the parameters in some range there are 
two possibilities for approximative procedures. 

3.1 The Iterative Solution of the Gap Equation 



Equation ( ^Tf ) can be iterated. For this we write it in the form M = ^JVhJ^M). 
Inserting some zeroth approximation M° into the rhs. the next approximation 
can be obtained which in turn must be inserted into the rhs. and so on. This 
procedure converges until the module of the derivative of the rhs. with respect to 
M does not exceed unity. This derivative takes its largest value at M = 0. From 
this a sufficient criterion for the convergence can be derived. Let us consider the 
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Figure 2: Plot of the rhs. and the Ihs. of the gap equation (pTj) for $c = 0.6, 
A = 1, T = 1 



derivative 



dr 



hs 



dM |„=o 
It is smaller than unit for 



-3AT/47r 



A < 



4m2 327r2$ 2 Stc^ 



3T2 



3T2 
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4m2 327r2$^2 



\ V3T2 



3T2 



9 



2^2 



Also, even for A obeying this restriction, the denominator may become small and 
make the iteration procedure diverge. In fact, this happens for M close to zero. 

In case the sufficient conditions are satisfied the iteration procedure works 
well. It does so even in the complex region, i.e., for M'^ < and down to T~0. 
The solution appears in powers of a/A. 



3.2 The Approximative Solution 

Instead of iterating the Eq. (^) one may approximate it and solve the approx- 
imated equation. We restrict ourselves to temperatures of order T~A~2. This 
is just the order where the phase transition happens (see below). Then in the 
gap equation the approximation ( pT] ) may be used. So the approximated gap 
equation reads 



T2 



-m^ + 3A$/ + 3A <^ — 



MT 



12 Atc 



+ 



167r2 



With the notation = 1 - |^ [in ^^^g^ - 27 



2m2 



the solution reads 



(33) 



M 



3AT ^-1 



3AT 



m25i + 3A<l>c^ + 



AT2 



(34) 
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Again, it can be seen that it is real for > (^^Si - 12$^^) / (l + 
Thereby, for < this is the case for > - T'^) whereas for 

> the solution is real for any <l>c. For + 3A<l>c^ + XT'^/A < the 

solution becomes complex and the mass squared becomes negative. The solution 
M^($c ) of the gap equation is shown in figure |^ for a temperature where it has 
an imaginary part for small $c and is real for larger $c- It is a nearly linear 
function and follows mostly the tree solution, M = fx, fx^ = —m? + 3A$c^- Only 
near M = it shows a nontrivial behavior which is the non perturbative region 
for small M discussed above. 

The validity of the approximation leading to Eq. (^) can be established 
as follows. First, consider large T. The solution is 

Second, consider T ~ l/VX. The solution is M^\/~\T. So in both cases the 
contributions neglected in (^) are of order ^ ~ or smaller (for example, at 
— m^(5i +3A$c^ + = the solution M = is exact). Hence including higher 
order terms into (0) would result in corrections of order V^. On the other hand 
side dropping any of the terms in (|6l|) changes the solution quantitatively. It can 
be checked that in doing so even the order of the phase transition changes. 




Figure 3: Real part (solid line) and imaginary part (broken line) of the effective 
mass squared as function of solving the gap equation for A = T = 1, the 
magnification shows details near M = 
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4 The Shape of the Effective Potential and its 
Extrema 



The effective mass and the effective potential can be calculated numerically and 
plotted. In the temperature region where the phase transition takes place the 
approximative solution (|3^ of the gap equation together with formula (|57|) makes 
it very easy to plot the effective potential as a function of $c, see figure 0. It shows 
the typical behavior of a first order phase transition. This important property 
can be shown exactly. 




Figure 4: The effective potential as function of $c for several temperatures in the 
region T_ < T < for A = 1 



We investigate the extrema of the effective potential. Differentiating V^s as 
given by Eq. (|30|) we obtain 



dVr 



eff 



+ 1 



6A 



+ + -Ao(M) 



1 



9$, 



(35) 



where (^Tj) was used. By means of the gap equation ( P?]) this expression simplifies 
and from the condition dVes/d^c^ = we obtain 



A$, 



(36) 



This relation must be considered together with the gap equation where we elim- 
inate A$^ by means of Eq. ( JI7\ ) and obtain 



= 2m^ - 6AAo(M) 



(37) 
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as the gap equation in the extrema of the effective potentiaL 

Eq. (|37D has a form similar to the gap equation (p7[) . Again, we consider the 

graphical solution written in the form M = \j2m? — 6AAo(M). Typical pictures 
are shown in figure |. 




Figure 5: Plot of the rhs. and of the Ihs. (= M) of the gap equation in the 
minimum as function of M for A = 1 and several temperatures 

The curve for the rhs starts in 2m? 5i — with a positive derivative (cf. 

Eq. (pT|)). Therefore, for sufficiently small T there is one real solution. This is the 
position of a minimum of the effective potential (for T = one obtains the known 
picture of the tree approximation). For > Am'^Si/X the starting point of the 
curve representing the rhs. is below the origin and there are two real solutions. 
Because the effective potential is a continuous function, this second solution must 
be a maximum. It appears at smaller values of $c than the minimum. So 

T ^ (38) 

is the lower spinodal temperature. Raising the temperature lowers the curve 
representing the rhs. further. The two solutions become closer to each other 
until they merge finally. At that temperature the minimum and the maximum 
are at the same place. Hence this is the upper spinodal temperature T+. For 
even higher temperature there is no real solution and the effective potential is a 
monotone function - the symmetry is restored. 

Because these are exact, non perturbative properties of Eq. 
exactly that the phase transition is of first order. 

It is important to note that the effective action is real in both its minima. 
For T < T_ there is only one minimum (that at finite $c)- There the effective 



37) it is shown 
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potential is real. This is clear because the effectve mass is real as shown just 
above and because the functions Ao(M) and Vi{M) are real for real arguments. 
In $c = the effective potentiual has a maximum and it has an imaginary part. 
This corresponds to the instability of this state. 

Consider T_ < T < T+. Here the effective potential has two minima, one in 
$c = and the other one at finite $c- In both minima (as, in general, for any 
T > T_) the effective potetential is real. So in the given approximation they 
are both stable. In fact, the effective potential should have an imaginary part 
in the minimum at $c = for T < Tc because there is a finite probabiity for a 
tunneling transition to the lower lying minimum. But this is beyond the given 
approximation. 

It is interesting to note that in T = T_ the effective mass at = vanishes. 
Also the second derivative is zero there, see Eq. (^if) . So this point is metastable 
as it must be in a first order phase transition. 

These features can be seen explicitely in the expansion of Eq. (^) for small 
A. Using ( pT| ) in the rhs. of Eq. ( |57[ ) we obtain 



= 2m^ - 6A 
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MT 1 
+ 



Att 167r2 



In 



2ni? 



+ 2m' 



(39) 



which can be easily resolved with respect to M. With the notations for 5i intro- 



duced in the preceding section and 83 = 1 + ^ {\a 27^ we obtain 



SAT 



Min/M ax 



± 



An 63 \ \An6: 



SAT 



+ 



2m^i AT2 



63 



263 



(40) 



where the upper sign corresponds to the minimum and the lower sign to the 
maximum of the effective potential. The condensate $c is related to this mass by 
Eq. (p6|). The lower spinodal temperature T_ is the temperature at which the 
maximum appears at Mm^x = when raising the temperature. It is the same 
as given by formula (pSf). Also, the upper spinodal temperature T+ follows from 
formula (^Of). It is the temperature for which the positions of the maximum and 
of the minimum coincide, i.e., where the square root vanishes: 



2m 



Si 



9A 

87r2 53 



(41) 



Within the given precision^, i.e., in the lowest nontrivial order in A, these tem- 
peratures can be written as 



T_ 



2m 

71 



SA 
16^ 



T^ 



SA 



2m I 

71 V ^^IStt^ 



(42) 



•^bearing in mind that corrections of order and higher foUow from the graphs contained 
in W^iPi(O), Eq. (|29|) 
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They are both of order m/\/~\ justifying the approximation made in the gap 
equation. The interval between T_ and T+ is of order ^fXm. 



X=0.1 




0.4 - 



0.2 0.4 0.6 0.8 1 vAT/2 

Figure 6: = M/V2 in the minimum (resp. maximum) (upper resp. lower 

part of the curves) of the effective potential as function of the for several 

values of the couphng A 

The critical temperature T^, which is the temperature where the effective 
potential in its minimum at finite $c equals its value at $c = 0, is in between 
these two temperatures, closer to T+. 

The positions of the minimum and the maximum of the effective potential as 
function of the temperature can be calculated from Eqs. (0) and ( PBD numeri- 
cally. Also, they can be taken from (|40|) except for small values of T where the 
approximation made in (|39|) is not valid (there the solution Mof Eq. (^) has a 
nonanalytic behaviour of the type M^'^^^/2m + 0(exp(— ^^^^))). The result is 
shown in figure ^ for several values of the coupling A. 

It is possible to calculate in a quite easy way the depth of the minimum at 
T = T_. For this we need the potential at $c = for this temperature. There 
the gap equation has the (exact) solution M = 0. This is just the situation when 
the curve representing the rhs. in figure § starts from zero. Than by means of 
and (|2^) we obtain 

Vesi^^^o, T^T_ 90 ■ 

In order to calculate the effective potential in its minimum at T = T_ we use 
OD and insert $c from Eq. (^). We obtain 
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Now we insert T_ from Eq. ( ^21) and MMin from (^). Expanding Vi according 
to (|57D we obtain for the difference 



VcS\Min, T=T. - Kffls^^o. T=T_ - ~~[Q~K^ ' ''^'^^ 

where higher corrections in A are dropped. This is the depth of the minimum at 
the lower spinodal temperature. 

It is interesting to calculate the effective mass M in the minima of the ef- 
fective potential as function of the temperature. This is the so called Debye or 
temperature dependent mass. For <l>c = it follows from Eq. (|3^). At T = T+ we 
note the value Mq = ^{■\/?> — l)\/~\m. In the other minimum it is what we denoted 
by MMin, given by Eq. (^) (upper sign). At T = we note MMin = 2m^ which 
is the tree value. At T = T_ we have MMin = 3-\/Am/27r and at T = we have 
MMin = 3v^m/47r. Both are shown in figure |^ as function of the temperature. 



M . 

mm 




T T 

+ 



Figure 7: The Debye mass, i.e., the effective mass in the positions of the minima 
of the effective potential as function of the temperature for A = 1 



Another important check of the consistency of the given approximation is 
that in the minimum of the effective potential the term linear in r] in the effective 



action (p8|) vanishes exactly. This can be seen by inserting $c from (|36| ) and 
using Eq. (^7|). Note that from this the disappearance of all tadpole diagrams in 
the Green functions, not only in the effective action, follows. 

Also, the connection between the effective mass and the second derivative 
of the effective potential in the minimum can be checked. For this reason we 
calculate c}^\4ff/(?(^c^)^ from (|35|) and obtain by means of (^) and (|37|) 



g^Kff _ 1 + 3AA^(M) 
d{^,f ~ 1 - 3AA'o(M) 



M^ , (44) 



with A' = {dA/dM'^) which is the expected relation within the given approxi- 
mation. 
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5 Comparison with different approaches 



Since the pioneering papers back to the 70th there have been several steps in the 
developement of the topic. Aheady in |^ the necessity of summing 'daisy' and 
'super daisy' diagrams was underhned. Perhaps the first summation of 'daisy' 
diagrams (also called 'ring' diagrams) for different models of quantum field the- 
ories was carried out in |T^, later on it was reconsidered in |]18| with specific 
application to the electroweak phase transition in the Standard Model. In the 
latter paper using some high temperature approximations a first order phase 
transition was obtained where also the imaginary part comes out right. As we 
will see below, this is a qualitatively good approximation. Detailed information 
about properties of the high temperature daisy approximation can also be found 
in IT 



5.1 One Loop Approximation 

The simplest approximation is the 'pure Tr In', which means the one loop con- 
tribution to the effective action without summing up any 'daisy' diagrams. The 
corresponding expression is simply 

Kff = + ^*c^ + Vl(/i^) (45) 
with = — + 3A$c^ and Vi{m) = |Tr ln(p^ + m^). The function Vi is given 



in the Appendix, (0), (0). The behaviour of this effective potential is in general 
well known, but there are some subtles. In order to encounter them we consider 
the extrema of (^) by equating zero its derivative 

^ = -^+2^c^ + 2AAo(/.^)> (46) 
whereEq. (^) had been used. This can be rewritten in the form 

/i^ = 2m^ - 9AAo(;u) . (47) 



This relation is like Eq. (|37D with the difference that the effective mass is /i^ = 
— -|- 3A$c^ in this case. In complete analogy to Eq. (^) the solution can be 
investigated graphically. Figure ^ applies changing merely the coupling according 
to A ^ |A. Hence the effective potential has a minimum for < 8m^/3A. For 
> 8m^/3A there are both, a minimum and a maximum. Raising T further 
they merge and disappear. These are properties of real solutions, i.e., for /i^ > 0, 
i.e., A$c^ > rn? 1^- For A$c^ < 'm.^/3 the effective potential is complex. Therefore, 
there might be another minimum if considering only the real part. The effective 
potential resp. its real and imaginary parts can be easily plotted using formula 
and the integral representation in (|55|), see figure ^. 
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Figure 8: The effective potential in one loop approximation for A = 1. Curve 
la resp. lb shows the real resp. imaginary part for a temperature where the 
effective potential is real in the lower minimum, curves 2a resp. 2b are for a 
higher temperature and the effective potential is complex in the lower minimum. 




Figure 9: The effective potential in one loop approximation for A = 0.1. Curves 1 
resp. 2 show the real resp. imaginary part, the magnification shows both minima 
present. 
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The general behaviour is that of a second order phase transition. In a narrow 
temperature region ~8m^/3A the real part has two minima. At < 8m^/3A 
the left one is lower and for > 8m^/3A the right one is the lower one. In 
the left one the effective potential is complex and in the right one it is real. It 
seems that the very tiny dip which is shown in figure § for A = 1 and in figure 
1^ for A = 0.1 in the magnification had not been observed before. It is clear that 
the pure one loop effective potential is unphysical because it is complex in the 
minimum at some temperatures. The position if the extrema are shown in figure 



ll| , curves 3,4 and 5 as functions of T. Curve 3 shows the minimum where is 



real, on curve 5 it is complex and curve 4 shows the position of the maximum. 

5.2 'Simple Daisy' Diagrams Summed Up 

After summing up all 'daisy' and 'super daisy' diagrams, which resulted in the 
effective potential (|29|) and the gap equation (^71) it is quite simple to go the 
step back to the simple 'daisy' diagrams summed up by substituting the tree 
value fj? = —m^ + 3A$c^ instead of the effective mass at two places in these 
formulae which results in 

Kff = -^^c' + 4$c^ - ^AAo(/i) + Vr{M) (48) 

with 

= -m^ + 3A<l>c^ + 3AAo(/i) , (49) 

where Vi resp. Aq are the same functions ([53| ) resp. ( pS] ) as in section]^. However, 
we have to show that (^8]) and (|^) is just the result of summing the simple 'daisy' 
diagrams. For this we have to go back to formula (^ and must pick up the daisy 
diagrams from 1^^^^. We obtain 



^ _!!!!$ 2^ A 4 1^^ in(p2+^2)_l C)_l p_l . . . (50) 

2 4 2 ^ 8 O 16 O 48 OO ^ ' 

To be explicit, here we used the notation Tr = Y4L-00 I Tr ln(p^ + 

m^) = — Tr A = — Vi(m) where A = l/(j9^ + m^) is the line in the graphs 
(in momentum representation). The vertex factors are —6 A. Now, the analytic 
expression for the graphs reads 



vb ^ _!!!!$ 2 A 4 1 

eff 2 4 2 



ln(p2 + m^) (51) 

1 -6AAo(m) _ 1 (-6AAo(m))^ _ 1_ {-QXAo{m)f _ 
"4 p2 + m2 8 (p2 + m2)2 24 (p2 + ^2)3 
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where Ao(m) = Tr + m^) is in fact the same as (^) resp. (|T^) written for 



T 7^ and with m instead of M. The graph 'eight' 
two ways 



-6A Tr 



-6A Tr 



S , can be represented in 
Ao(m) 



p2 _|_ ^2 



(cf. Eq. (|T2D). Eq. ( pT)) is aheady the expansion of the logarithm (ln(l + e) = 
e — + — • • •) except for the 'eight', which enters with the wrong sign. In 
writing the coefficient in front as | = — | + 1 (times 3 A after cancelhng a factor 
of 2) we obtain 



cff 



m ^ 



2 



4 



+- Tr 
2 



ln(p^ + 



1 3AAnM 



+ In 1 + 



3AAo(m) 

p2 _|_ ^2 



2-$,2 + -<|.,^--AA^(m 



2 + 

-Tr ln(/ + m2 + 3AA, 



o(m 



This is just 



and (M). 




Figure 10: The effective potential with 'simple daisy' graphs summed up for 
A = 1. Curves la, 2a, 3a resp. lb, 2b, 3b show the real resp. imaginary parts for 
T = 1, 1,4, 3.4. 

The investigation of the extrema of is more complicated that in the other 
cases and not very interesting because this case in anyway not physical. Never- 
theless it is interesting to comment on them. can be plotted qute easily using 
the above formulae. It is shown in figure |1^ as function of for several values 
of the temperature. It has an imaginary part for A$c^ < ?Ti^/3. At A$c^ = rn^/'i 
it has a cusp and above it is real. Again, the real part has two minima in a 
narrow temperature range. In general, the behaviour is quite similar to that of 
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Figure 11: Positions of minima and maxima of the effective potential at A = 1 in 
different approximations, see text 

tlie 'simple Tr In'. The positions of the minima can be found numerically and 
they are shown in figure |ll] as curves 6 and 7 (the maximum is not shown there). 
In the minimum represented by curve 7 resp. 6 the effective potential is complex 
resp. real. 

5.3 High Temperature Approximation and higher loop 

Frequently high temperature approximations for the functions Vi and Aq are used 
which are given by parts of the Eqs. (|57|) , (|6T|) . Another type of approximations 
used can be obtained by taking only the / = contribution in the sum over the 



Matsubara frequencies. A typical example is given in |jT9[. There, the expression 



= + \x^c' +'^i-m' + 3A$.^) -Y^[^-m' + 3A$e^ (52) 



has been obtained (see also fT^]). Here, the phase transition is first order, the 
effective potential is real in its minima. The positions of the maximum are shown 
in figure |TT] ( curve 8 resp. curve 9). It is seen, the behaviour is similar to that 
of the correct solution (curves 1 and 2), but the numbers are about of 10 ... 20 
% higher. So we conclude that the simple formula (^) works basically well and 
gives a result which is qualitatively correct and quite close to the right numbers. 



In the paper the daisy and super daisy graphs have been summed up by 
direct summation of the graphs performing the difficult task of calculating the 
combinatorical coefficients (which in the present paper was done implicitely in 
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section A first order transition resulted. Than the next graph, , had 

been included in the resummation. However, the combinatorics could not been 
done to the end. This is probably the reason that an unwanted imaginary part 
appeared which led the authors to the conclusion that the subleading terms they 
calculated are not consistent. 



In another paper a similar approach had been undertaken trying to solve 
the gap equation with two graphs. In a subsequent unpublished paper also, 
a first order transition was proposed. 



6 Conclusions 

In the foregoing sections the effective potential for a scalar theory with sponta- 
neous symmetry breaking and finite temperature is calculated. All 'daisy' and 
'super daisy' diagrams are summed up and the properties of the resulting gap 
equation are investigated. It is shown exactly that the phase transition is first 
order. The imaginary part of the effective action comes out in the correct way. 

It is important to notice that higher loop corrections cannot change these 
features. The point is that they are all expressible in terms of Feynman graphs 
with propagators l/(/c^ + M^) (in momentum space) where the mass M is solution 
of the gap equation. As shown it is positive in the physical region, in the minima 
of the effective potential for instance (see figure |^. Consequently all graphs 
have lines with non vanishing denominator. As the theory is Euclidean, no more 
infrared divergencies can occur. A special discussion deserves the contribution of 
the graph . As the vertex factor is A$c and in the minimum of the effective 

potential we have \/A$c ~ 1 it seems to contribute the same order as the daisy 
graph. However, this cannot change the conclusion drawn above. It is only the 
quantitative characteristics of the phase transition which may become altered, 
not the qualitative ones like its order. 

Natural extensions of the present work is the inclusion of fermions, gauge 
bosons and external fields. 
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Appendix 

The function Vi(M) appearing in (p9D is the notation for the 'Tr In' term: 

Vi(M) = --Tr lnA(M) = -^5: J J^s {^f + + M') (53) 

with (3 = 1/kT and uJi = 2ttI/i3. Removing the ultraviolet divergencies and taking 
into account the normalization conditions the explicit expression reads 



ViiM) 



647r2 



In 



2m? 



M^T' (M\ 



27r2 ^\Tl 



(54) 



where the last term is the temperature dependent contribution. The function 5*2 
can be represented as a fast converging sum over the modified Bessel function 
K2 and by an integral representation which is suited for complex values of the 
argument as well: 



°° 1 1 /■ 

-^2(3;) = Y. — K2{nx) = — / 



(^2 _ ^2)3/2 

an 

>x e" - 1 



Its expansion for small x reads 
^4 



^ , , TT' TT^ TTX X^ ( 3 , X \ ^, q, 

S^ix) = + — + — 27-- + 2n— + 0(x^^ 

^ ' 45a;2 12 6 32 V 2 Att } ^ 



(55) 



(56) 



giving rise to the high temperature expansion of Vi which is at once the expansion 
for small M: 



Vi{M) 



+ 



90 
1 



24 



647r2 



\2 



(57) 



1 (4vrT) 
In ^ ^ - 27 



2rTi2 



Here 7 is the Euler constant. 

Similar formulae hold for the 'daisy' graph. Although they can be derived by 
means of (|3l|) it is useful to have them separately: 



Ao(M) 



I ^ r ^ 1 



167r2 



2nf + 



In 



2m? 



1 



where, again, the last term is the temperature dependent contribution with 

00 1 1 



5*1(0;) = ^ — Kiinx 



n=l 



n 



dn 

X Jx e 



— x^ 



(59) 
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The corresponding expansions are 
and 

^ S - ^ - lib {-1- ^ -'"1 — (f ) ■ 

(61) 
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